Exploring binary-neutron-star-merger scenario of short-gamma-ray bursts by 

gravitational-wave observation 



O 

(N 

m 



43 : 

Oh- 
i < 

o ■ 

c3 . 



> 
ON 
OO 

(N 

(N 
O 
O 



Kenta Kiuchi, 1 Yuichiro Sekiguchi, 2 Masaru Shibata, 3 and Keisuke Taniguchi 4 

'Department of Physics, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan 

Division of Theoretical Astronomy /Center for Computational Astrophysics, 
National Astronomical Observatory of Japan, 2-21-1, Osama, Mitaka, Tokyo, 181-8588, Japan 
3 Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto, 606-8502, Japan 
^Department of Physics, University of Wisconsin-Milwaukee, P.O. Box 413, Milwaukee, Wisconsin 53201 

(Dated: February 13, 2010) 

We elucidate the feature of gravitational waves (GWs) from binary neutron star merger collapsing 
to a black hole by general relativistic simulation. We show that GW spectrum imprints the coa- 
lescence dynamics, formation process of disk, equation of state for neutron stars, total masses, and 
mass ratio. A formation mechanism of the central engine of short 7-ray bursts, which are likely to be 
composed of a black hole and surrounding disk, therefore could be constrained by GW observation. 
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Introduction. — Coalescence of binary neutron stars 
(BNSs) is one of the most promising sources for 
kilometer-size laser interferometric gravitational-wave 
detectors It is also a candidate for the central en- 
gine of short 7-ray bursts (GRBs), which emit huge en- 
ergy > 10 48 ergs in a short time scale ~ 0.1-1 s [1, 
According to a standard scenario of GRBs based on the 
so-called merger scenario, a stellar-mass black hole (BH) 
surrounded by a hot and massive disk (or torus) should 
be formed after the merger. Possible relevant processes 
to extract the energy of this BH-accretion disk system 
for launching a relativistic jet are neutrino-anti neutrino 
{vv) annihilation [J] and/or magnetically driven mech- 
anisms, so-called Blandford-Znajek process [5(. Studies 
based on the vv annihilation scenario suggest that an ac- 
cretion rate of M > O.lM /s is required to achieve a suf- 
ficiently high energy efficiency @ and that if the disk had 
a mass > 0.01M©, it could supply the required energy by 
neutrino radiation for duration of > 100 ms [7|. Gen- 
eral relativistic magnetohydrodynamic simulations also 
indicate that if a rapidly rotating BH is formed, its rota- 
tional energy can be extracted by the Blandford-Znajek 
process to achieve a high energy efficiency 0. However, 
the merger scenario has not been proven yet observation- 
ally, because the counter part for most of short GRB has 
not been identified [9|. Gravitational waves (GWs) are 
much more transparent than electromagnetic waves with 
respect to scattering with matter, and hence, they can 
propagate from extremely dense region with negligible 
scattering. Therefore, GWs can be valuable observable 
for determining the central engine of GRBs [l(| and for 
exploring its formation process. This fact motivates to 
study coalescence of BNSs and emitted GWs. 

For theoretically studying the late inspiral, merger, 
and ringdown phases of BNSs, numerical relativity is the 
unique approach. Simulations of BNS merger 11- 13j| 
suggest that (i) if the total mass of the BNS is smaller 
than a threshold mass (-Mthr) a hypermassive neutron 



star (HMNS) sustained by rapid and differential rotation 
is formed after the merger and will survive for more than 
~ 100 ms [Hj], and (ii) otherwise, a BH is formed after 
the merger in a dynamical time scale ~ 1 ms. In the for- 
mer case, GW emission [ll[ or magnetic-field effect such 
as the magnetorotational instability [3| or neutrino ra- 
diation [1 51 ] will play an essential role in the subsequent 
evolution of the HMNS, which will eventually collapse 
to a BH. Because of the complexity, it will not be an 
easy task to theoretically clarify the formation process of 
the BH and gravitational waveforms as well as to accu- 
rately predict the launching mechanism of short GRBs. 
By contrast, in the later case, the formation process of a 
BH and surrounding disk is primarily determined by the 
general relativistic hydrodynamics, and magnetic field or 
neutrino radiation will play a minor role. Reliable the- 
oretical prediction of merger process and gravitational 
waveforms by numerical relativity is feasible in this case. 
Motivated by this idea, we systematically performed nu- 
merical relativity simulations to clarify the relation be- 
tween possible formation process of GRB's central engine 
(a system composed of a BH and surrounding disk) and 
gravitational waveforms, and suggest that GW observa- 
tion could constrain the formation process for the central 
engine of short GRBs. 

Method and initial models. — Formulation and numer- 
ical schemes for solving Einstein's equations, hydrody- 
namics, and other techniques such as GW extraction 
are essentially the same as those in Ref. [l3|, to which 
the reader may refer for details. For modeling the EOS 
of neutron stars, we adopt the Akmal-Pandharipande- 
Ravenhall (APR) 16], Skyrme-Lyon (SLy) [l^, and 
Friedman-Pandharipande-Skyrme (FPS) [18| EOSs for 
zero-temperature part. Shock-heating effect, although 
it gives only minor contribution, is mimicked by T-law 
EOS with T = 2. BNSs of the irrotational velocity field 
and zero-temperature in quasiequilibrium circular orbits 
are prepared as initial conditions, which are computed 
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FIG. 1: + modes of GWs for (a) A2.9-0.8, (b) S2.8-1, and (c) 
F2.6-0.8. D is the distance from the source to the observer, 
who is located along the axis perpendicular to the orbital 
plane. Emerge denotes approximate time at the onset of the 
merger. 



by the LORENE library codes [19(. For all the mod- 
els, the initial orbital separation is large enough that the 
BNSs spend about 4 orbits in the inspiral phase before 
the merger. Grid structure and resolution are approxi- 
mately the same as in Ref. [131 ]; the major diameter of 
the neutron stars is covered by about 80 uniform grids 
and outer boundaries of numerical domain are located at 
about 1.2Ao where Ao is initial gravitational wavelength. 
Total mass, too, is chosen to be large enough that a BH 
is formed in a dynamical time scale (~ 1 ms) after the 
onset of the merger. This can be achieved for too > M t hr, 
where M thl « 2.8-2.9M Q for the APR, « 2.7-2.8M Q for 
the SLy, and « 2.5-2.6M Q for the FPS EOS, respec- 
tively [ill ]. Total mass and mass ratio are chosen in the 
range of 2.9M Q < to < 3.1M ,O.8 < v = m 1 /m 2 < 1 
for the APR, 2.8M Q < m < 3.0M o , v = 1 for the SLy, 
and 2.6M < to < 2.8M Q ,v = 0.8 and 1 for the FPS 
EOS. mi and TO2 (to.2 > toi) denote the gravitational 
mass of two neutron stars in isolation and mo = mx+m?,. 
We name the models according to the EOS, mo, and 
v; e.g., A2.9-0.9 is modeled by the APR EOS with 
too = 2.9Af Q and v = 0.9. 

Results. — We investigate the relation between GWs 
and possible formation process of a BH-disk system in 
a wide range of parameter space. First of all, we briefly 
review a typical coalescence process of BNSs collapsing 
dynamically to a BH (see Ref. 13j in details) . During the 
inspiral phase, the BNSs adiabatically evolve gradually 
decreasing its orbital separation due to the gravitational 
radiation reaction. GWs in this phase are characterized 
by a chirp signal (see the waveform for i ret — imergo 
in Fig. [IJ. After the binary separation becomes < 3 
neutron star radii, the merger sets in. For v ~ 1, a 
high-density region with density > 10 15 g/cm 3 is sub- 
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FIG. 2: Effective amplitude for models A2.9-1, A2. 9-0.8, S2.8- 
1, F2.6-1, and F2.6-0.8. The filled circles denote the QNM fre- 
quency of the formed BHs. The small panel shows a schematic 
figure of the GW spectrum. 



sequently formed at the center, and then, collapses to a 
BH within ~ 1 ms. During this merger process, small spi- 
ral arms, composed of two symmetric arms, are formed 
around the BH but they are eventually swallowed by 
the BH. The final outcome is a rotating BH nearly in 
a vacuum state. For v < 0.9, the less massive star is 
tidally deformed in a close orbit and disrupted just after 
the onset of the merger. The outer part of the tidally 
disrupted neutron-star matter forms an asymmetric spi- 
ral arm which is more massive and widely-spread than 
that for the equal-mass case. Although the central high- 
density region, formed at the onset of the merger, col- 
lapses to a BH, the large spiral arm subsequently forms 
an accretion disk around the BH because it has angular 
momentum large enough to escape from the capture by 
the BH. GWs in the merger phase are characterized by 
the short-term burst-type waves (0 < t rct —t mclgc < 1 ms 
in Fig.QJ and after the BH formation, quasinormal modes 
(QNM) are excited (t rc t — t mclgc > 1 ms in Fig. [T]) . Fig- 
ure[fl which plots GWs for A2.9-0.8, S2.8-1, and F2.6-0.8 
and shows typical waveforms for the BNS merger, illus- 
trates features mentioned above. 



GW spectrum more clearly reflects dynamics of the 
merger process. Figure [2] plots an effective amplitude of 
GWs as a function of the frequency, /i e ff(/), for models 
A2.9-0.8, A2.9-1, S2.8-1, F2.6-1, and F2.6-0.8. The effec- 
tive amplitude is defined by h e s = h(f)fmo/D with h{f) 
being the Fourier transform of h+ + ih x . We assume that 
the distance from the observer to the source is D = 100 
Mpc, because more than one event per year is predicted 
for such a large distance [2(j. Note that GW frequency 
increases as the orbital separation decreases and as the 
compactness of the merged object increases: GWs with 
/ < 1 kHz, 1 < / < 3 kHz, and / > 3 kHz are emitted in 
the inspiral phase, early merger, and late merger phases, 
respectively. 
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FIG. 3: G/cutf^i/c 3 as a function of compactness of less mas- FIG. 4: Disk mass as a function of CT"// pea k for all the models, 
sive binary component for all the models. Sequence A(l) rep- 
resents the models A2.9-0.8, A3-0.8, and A3. 1-0.8 from left to 
right. In a similar way, A(2) is a sequence of A2.9-0.9, A3-0.9, 
and A3.1-0.9, A(3) is A2.9-1, A3-1, and A3. 1-1, S(l) is S2.8-1, 
S2.9-1, and S3-1, and F(l) is F2.6-1, F2.7-1, and F2.8-1. The 
models below the horizontal line of G/ cu tfni /c 3 = 0.02 can 
produce disks of mass > 0.01M©. 



The effective amplitude gradually decreases for / < 3 
kHz according to oc / _ " with n w 1/6 which is the typical 
number for the inspiral phase, and for / > / cut ~ 2.5-4 
kHz, its gradient becomes steep (see the small panel in 
Fig. [2] for / > / cut ). For / < / cut , the merged object 
has a binary-like structure (i.e., there exist two density 
maxima) which enhances emissivity of GWs, whereas at 
/ ~ /cut, such structure is disrupted, resulting in a quick 
decrease of the GW amplitude. A hump, which appears 
for 4 < / < 7 kHz with the peak frequency / pca k ~ 5-6 
kHz, reflects the formation and evolution of the spiral 
arms orbiting the central object. For / > / pca k, ^cff 
exponentially decreases; this is the typical feature for 
the spectrum associated with a QNM ring-down of the 
formed BH. These spectrum features are qualitatively 
universal as shown in Fig. [2J However, the spectrum 
shape depends quantitatively on the EOS, mo, and v. 

The spectrum shape for / > / cut is characterized 
by /cut, /peak, and peak amplitude and width of the 
hump (/ipeak and a), for which a schematic figure is 
described in the small panel of Fig. [2] We perform a 
fitting procedure to extract these characteristic quanti- 
ties from the spectrum. As the first step, the spectrum 
is divided into two parts; one is for / < fa, inspiral 
and early merger phases, and the other is for f > fa, 
the merger and ringdown phases, /d is typically cho- 
sen as 4 kHz. For / < fa, we fit the spectrum by 
/- 1/ %/[l + exp{(/-/ cut )/A/}]. For / > fa, we adopt 
the Gaussian distribution /i pea kex 
the fitting, / cut , A/, h , h 
mined by the x-square fitting [2l|. 

Figures [3] plots G/ cu t"T-i/c 3 as a function of compact- 
ness of less massive companion defined by Gmi/c 2 R 



where R is its circumferential radius. G/ cu tmi/c 3 corre- 
lates strongly with Gmi/c 2 R. G/ cut mi/c 3 depends also 
on mo; for the larger value of mo, it is larger for a given 
EOS and v. Note that larger values of mo and v result in 
smaller values of disk mass (see Fig. 4). Thus, / cu t is an 
indicator that shows the occurrence of the tidal disrup- 
tion and possible disk formation. (We note that to obtain 
/cutfii, we have to determine mi from the gravitational- 
wave signal in the inspiral phase |22j . It is also worthy to 
note that R may be determined, if / cu t is measured and 
mi is determined; see Ref. 23j for the related topic.) 

Figures [4] plots the mass of disk surrounding a BH as a 
function of cr// poa k. This shows that the disk mass has a 
positive correlation with a/ /peak- On the other hand, we 
find that the disk mass does not have a clear correlation 
with /ipeak- The disk mass also correlates with mo — M t hr 
and v. mo — Mthi < 0.2M Q and v < 0.8 appear to be 
necessary for producing a massive disk with > O.OIM© 
(see A2.9-0.8, F2.6-0.8, and F2. 7-0.8). The disk mass for 
other models which do not satisfy this condition is 10~ 3 
5 M m " 
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k exp[-(/-/ P cak) 2 / cr2 ]- In 
^eakj/peak, and a are deter- 



[l3j ]. Combination of Figs. 3 and 4 also 
proposes that a necessary condition for producing a disk 
of mass > O.OIMq is G/ cu t"^i/c 3 < 0.02 and cr// pcak > 
0.25. 

The frequency of GWs in the merger phase is slightly 
outside the most sensitive frequency band of the first and 
second generation detectors such as LIGO and advanced 
LIGO (10-1 kHz), because the values of / cu t and / pca k 
are ~ 2.5-7 kHz. Nevertheless, the amplitude at / = / cut 
and /peak is not extremely small for a hypothetical dis- 
tance D = 100 Mpc. If the proposed third generation de- 
tectors such as Einstein telescope become available with a 
special design for high-frequency sensitivity, these quan- 
tities may be measured. Our results indicate that the 
formation and evolution processes of a BH and surround- 
ing disk are imprinted in GW spectrum. Also, GWs are 
the unique observable which directly carries information 
of the central region of GRBs. Thus, it is worthy to ex- 
plore methods for extracting physical information on the 
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merger of BNSs collapsing to a BH from observed GWs, 
which may clarify the formation process of short GRBs. 

Constraining the merger process by GW observation. — 
As mentioned in Introduction, the merger hypothesis for 
the central engine of short GRBs requires that the mass 
of disk around a BH is greater than ~ 0.01A/ Q . Our nu- 
merical results show that disk mass correlates strongly 
with cr//peak and G/ cu tmi/c 3 irrespective of EOS. As- 
suming that GWs from BNSs are frequently observed, 
and also, GWs and short GRB are observed simulta- 
neously for several events in the future, we propose a 
strategy for exploring the merger hypothesis by analyz- 
ing GWs. As a first step, we should determine mo, v, 
/cut, /peak, and a for each event: The values of mo and v 
will be determined from the inspiral waveform using the 
matched filtering technique [22| . and / cu t) /peak, and a be 
from the merger waveform. As a second step, we should 
infer the disk mass from G/ cu tmi/c 3 and cr// pca k. If 
the conditions of G/cut^i/c 3 < 0.02 and cr// poa k ^ 0.25 
are always satisfied for events in which GWs and short 
GRB are simultaneously observed, the merger hypothe- 
sis will be strongly supported. Because a small value of 
G/cut^i/c 3 (for the typical mass of neutron stars 1.3- 
1.5Mq) implies that mass ratio is much smaller than 
unity, we could conclude that the origin of short GRBs 
would be the merger of unequal-mass BNSs. By con- 
trast, if simultaneous detections occur irrespective of the 
values of G/ cu tmi/c 3 and cr// pe ak, this does not agree 
with our numerical results, because these should corre- 
late with the disk mass. In such case, we have to conclude 
that an unknown mechanism plays a role for producing 
the GRBs. 

Our numerical results indicate that massive disk is not 
formed for v w 1 or mo > 3M Q for any EOS. This im- 
plies that the merger of BNS collapsing to a BH does 
not always produce GRBs. For such case, only GWs will 
be observed. This suggests another test for the merger 
scenario. 

Finally, a comment is given for the case that an HMNS 
is formed. Because the HMNS should eventually collapse 
to a BH via gravitational radiation [11| or angular mo- 
mentum transport by magnetohydrodynamic effect 
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or neutrino cooling, the HMNS formation scenario can be 
an alternative in the merger hypothesis. In the HMNS 
formation, the GW spectrum for the merger phase is sig- 
nificantly different from that in the BH formation be- 
cause quasiperiodic GWs emitted by a quasiperiodic ro- 
tation of the HMNS is likely to produce sharp peaks in 
the spectrum and characterize the spectrum [1 11 ] . In this 
case, detecting these peaks of GWs will play a crucial 
role for confirming the merger hypothesis. 
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